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1. Introduction 

The origin of mass is the most urgent issue of the particle physics today. One of the candidates 
for the theory beyond the Standard Model towards that problem is the walking technicolor which 
is the strongly coupled gauge theory having a large anomalous dimension y m ~ 1 and approximate 
scale invariance due to the almost non-running (walking) coupling [jl|, g]. The walking behavior is 
in fact realized in the QCD with large number of (massless) flavors Nf which possesses Caswell- 
Banks-Zaks infrared fixed point (IRFP) ^ for the value larger than Nf ~ 8.0 in the two-loop beta 
function. The exact IRFP would be washed out by the dynamical generation of a quark mass m 
in the very infrared region \i < m for Nf < NJ, NJ being the critical number. However, for Nf 
very close to NJ, m could be much smaller than the intrinsic scale A (3> m), an analogue of Aqcd> 
beyond which the coupling runs as the asymptotically free theory, so that the coupling remains 
almost non-running for the wide infrared region m < fx < A as a remnant of the would-be IRFP. 
The case N > Nf" is called conformal window, although conformality is broken in the ultraviolet 
asymptotically free region beyond A. The critical number Nf" was estimated as NJ ~ 11.9 [|j] 
by comparing the two-loop IRFP value with the critical coupling of the ladder Schwinger-Dyson 
equation analysis []5j|. 

Although the above results from the two-loop and ladder approximation are very sugges- 
tive, the relevant dynamics is obviously of non-perturbative nature, we would need fully non- 
perturbative studies. Among others the lattice simulations developed in the lattice QCD would 
be the most powerful tool for that purpose. Our group, LatKMI Collaboration, was organized for 
such studies on walking technicolor as a candidate for the theory beyond the Standard Model. The 
immediate issues are: What is the critical number NjH What is the signatures of the walking theory 
expected to be slightly smaller than NJ1 The above two-loop and ladder studies suggest that the 
walking theory if existed would be in between Nf = 8 and Nf =12. The Nf = 8 in particular is in- 
teresting from the model-building point of view: The typical technicolor model [^] is the so-called 
one-family model (Farhi-Susskind model) which has a one-family of the colored techni-fermions 
(techni-quarks) and the uncolored one (techni-leptons) corresponding to the each family of the SM 
quarks and leptons. It can embed the technicolor gauge and the gauged three generations of the SM 
fermions into a single gauge group (Extended Technicolor) and thus is the most straightforward 
way to accommodate the techni-fermions and the SM fermions into a simple scheme. Thus if the 
Nf = 8 turns out to be a walking theory, it would be a great message for the phenomenology to be 
tested by the on-going LHC. 

To date, some groups [[7],[8],|9, 10, 11, |12|] carried out lattice studies on 8 flavors; Ref. [f7|] com- 



puted the running coupling constant in Nf = 8 by the Schrodinger functional method in staggered 
fermion case and reported that Nf = 8 is in the chiral broken phase. Refs. [Q, |9[ [K], [II , 12] inves- 



tigated the hadron spectrum with the standard Wilson fermions [Bp, the Stout improved staggered 



fermions JlQ|], the Asqtad improved staggered fermions Ql 1| ] and the naive staggered fermions [g, 
12]. The Ref. ^ concluded the Nf = 8 is in the conformal window, but other groups concluded 
that the Nf = 8 resides on the chiral broken phase. Even if Nf = 8 is in the chiral broken phase, it is 
not clear whether the behavior of this system is QCD-like or the walking with the large anomalous 
mass dimension. Nobody has investigated the possibility that Nf = 8 is in the walking. 

We simulate 8-flavor QCD with alternative lattice fermion, HISQ, in which the flavor symme- 
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try in the staggered fermion is improved and, of course, it is expected that the behavior towards the 
continuum limit is improved. We show the preliminary result of the hadron spectrum and analyzed 
the data based on the hyperscaling [13] as well as the chiral perturbation theory (ChPT). From 
the hyperscaling analysis, we derive the anomalous mass dimension y m . We observe qualitatively 
different behavior of the 8-flavor case in contrast to the 4-flavor case which shows clear indication 
of the hadronic phase as in the usual QCD. 



2. Simulation 



In our simulation, we use the tree level Symanzik gauge action and the highly improved stag- 
gered quark (HISQ) action without the tadpole improvement and the mass correction in the Naik 
term. See Ref. for the detail of the HISQ action. We use the MILC code [jl5|] with modifica- 
tions to simulate Nf = 4n HISQ by using the standard Hybrid Monte-Carlo (HMC) algorithm. We 
computed the hadron spectrum as the global survey in the parameter region and we obtained M K , 
M p , fjf and (yy) as the basic observable. 

The simulation for Nf = 4 is carried out at j8(= 6/g 2 )=3.5, 3.6, 3.7 and 3.8 for various quark 
masses on 12 3 x 16 and 16 3 x 24. We took over 1000 trajectories on the small lattice and about 
600 trajectories on the large lattice in Nf = 4 case. The simulation for Nf = 8 is carried out at 
P(= 6/g 2 )=3.6, 3.7 and 3.8 for various quark masses on 12 3 x 32 and 24 3 x 32, and for nif =0.02 
and 0.04 on 30 3 x 40 at j8 = 3.7. We took about 1000 trajectories on 12 3 x 32, between 200 and 
400 trajectories on 24 3 x 32 and 600 trajectories on 30 3 x 40. 



See Ref. [16] for our simulation in Nf = 12 and 16. 



3. Spectrum 

In this section, we show preliminary results of ChPT analysis and the finite-size hyperscaling 
analysis in 4- and 8-flavor cases. 

3.1 Nf = 4 

The result of Nf = 4 is shown in Fig. |l|. The data on 16 3 x 24 at j8 = 3.7, the pion mass 
squared, the decay constant and the chiral condensate are plotted on the panel from the left to the 
right respectively. M\ is proportional to rrif. f n and (yy) have the non zero value in the chiral 
limit. Thus, the Nf = 4 QCD has the property of the chiral broken phase and this is regarded as the 
signal of the chiral broken phase in the dynamical case of lattice QCD. 

3.2 Nf = 8 

In this subsection, we analyze Nf = 8 system by ChPT and the finite-size hyperscaling relation. 
In Fig. ||, as the typical example of the raw data, we plotted the preliminary data of M K and M p 
as a function of the quark mass ntf on 12 3 x 32 and 24 3 x 32. It is shown that M % and M p are 
similar behavior, that is, the plateau appears in ntf < 0.06 on small lattice (12 x 32) at all /3s and 
in nif < 0.02 on large lattice (24 x 32) at j8 = 3.8. In the following, we attempt ChPT and the 
finite-size hyperscaling analysis by using these data. 
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Figure 1: In Nf = 4 SU(3) gauge theory on 16 3 x 24 at /3 = 3.7; Left: M\ as functions of m/, Center: f K 
as functions of m/, Right: as functions of ntf. In all panels, the blue-solid line is the linear fit. 
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Figure 2: In Nf = 8 SU(3) gauge theory on 12 3 x 32 and 24 3 x 32 at J3 = 3.6, 3.7 and 3.8; Left panel: M n 
as functions of w/, Right panel: Mp as functions of w/. 



3.2.1 ChPT analysis 

We analyze the Nf = 8 data by ChPT. Here, we pick up the data at j8 = 3.7, which has the 
30 3 x 40 volume data. Fig. || shows M~, f n and (yy)- F° r ea °h the largest volume data is used 
for the quadratic-fit corresponding to ChPT analysis. The data corresponding to the plateau is not 
included into the fit. 

As the fit result, M\ in the chiral limit is consistent with zero. This is fairly independent of the 
fit range. The quadratic-fit of f % in the chiral limit indicates non-zero value. These properties are 
consistent with the chiral broken phase. On the other hand, the condensate is not inconsistent with 
zero in the chiral limit, in contrast to Nf = 4 case (compare Fig. |3| with Fig. [I]). 

The definite conclusion cannot be given here because x 2 /d.o.f. ~ O(10) in ChPT fit and 
because it is difficult to take the infinite volume limit. Even if Nf = 8 is in the broken phase, there 
is a possibility that the walking behavior as a remnant of the conformality can be observed in the 
form of hyperscaling relation. 

3.2.2 Finite-size hyperscaling analysis 

If the system is in the conformal window, physical quantities, M#, are described by the finite- 
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Figure 3: In 51/(3) gauge theory with 8-flavors on 12 3 x 32, 24 3 x 32 and 30 3 x 40 at J3 = 3.7; Left: M\ 
as functions of m/. % 2 /d.o.f. = 42. Center: f K . % 2 jd.o.f. = 105. Right: (ijryr). X 2 jd.o.f . = 154; In all 
panels, the solid line stands for the quadratic-fit. 



^{X) where JP(X) is unknown function with the scaling 



size hyperscaling relation | |13| ] ; LMh 

i 

variable X = Lml +Y . The 7 in this equation is defined as the anomalous mass dimension. We carry 
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out the hyperscaling analysis with our data of M# = {M n ,f n ,M p } by the following fit function 



LM, 



H 



1 

Co + cyLmy r 



(3.1) 



Figs. |5] and |6| are the finite-size hyperscaling result with this fitting of LM K , Lf n and LM p 
respectively. The filled symbol is applied to the hyperscaling fitting. These data on various lattice 
sizes align at an optimal value of 7. Although the fit quality is more or less the same level as that of 
the ChPTfit, we extracted the 7-value y(M n ) ~ 0.6, y{f K ) ~ 1.0 and y{M p ) ~ 0.8 at all jS, which are 
not universal. Therefore, our result of 8 flavors does not show the clear signature of the conformal 
window. Still, this situation in Nf = 8 is qualitatively different from the case of Nf = 4. Actually 
we did not find even an alignment of the hyperscaling for each physical quantity in Nf = 4. 






Figure 4: Left: The hyperscaling relation of LM K as functions of Lm\ ^inNf = 8 at j3 = 3.6, % 2 jd.o.f . 



41.7. Center: at J3 = 3.7, % 2 jd.o.f. = 39.8. Right: at j8 
The filled symbol is applied to the hyperscaling fit. 



: 3.8, X 2 jd.o.f. = 37.3. At all j8, y(M % ) - 0.6. 
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4. Summary 

We have made simulations of lattice QCD with 4 and 8 flavors by using the HISQ action. We 
obtained the following result; The Nf = 4 QCD is in good agreement with the chiral broken phase. 
The Nf = 8, on the other hand, does not seem to be inconsistent with both ChPT and the finite-size 
hyperscaling, with % 2 /d.o.f. not being small for both analyses. We extracted the /-value from the 
hyperscaling analysis, y{M n ) ~ 0.6, j{f%) ~ 1.0 and y{M p ) ~ 0.8 at all j8, which are not universal. 
Therefore, our result of 8 flavors does not show the clear signature of the conformal window. This 
may be an indication of the walking behavior that appears in the broken phase just below the edge 
of the conformal window. 

We should mention that there are several possible systematic uncertainties not considered in 
this report; Our data include those with slightly different aspect ratio, the ratio of the temporal 
length to the spatial length. Our result of the finite-size hyperscaling analysis may suffer from the 



systematic error due to the different aspect ratio. Furthermore, as pointed out in Ref. [17], there 



exists the mass correction in the hyperscaling relation for the heavy quark region. Also the linear 



ansatz adapted in Eq. (3.1) may not be sufficient to fit our data. To improve the situation for better 



understanding, we will accumulate more data for various fermion masses and j8's on larger lattices, 
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and carry out detailed analysis using those data. 
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